Transition to Landau Levels in Graphene Quantum Dots 
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We investigate the electronic eigenstates of graphene quantum dots of reahstic size (up to 80 nm 
diameter) in the presence of a perpendicular magnetic field B. Numerical tight-binding calculations 
and Coulomb-blockade measurements performed near the Dirac point exhibit the transition from the 
linear density of states at B = to the Landau level regime at high fields. Details of this transition 
sensitively depend on the underlying graphene lattice structure, bulk defects, and localization effects 
at the edges. Key to the understanding of the parametric evolution of the levels is the strength of 
the valley-symmetry breaking K—K' scattering. We show that the parametric variation of the level 
variance provides a quantitative measure for this scattering mechanism. We perform measurements 
of the parametric motion of Coulomb blockade peaks as a function of magnetic field and find good 
agreement. We demonstrate that the magnetic-field dependence of graphene energy levels may serve 
as a sensitive indicator for the properties of graphene quantum dots and, in further consequence, 
for the validity of the Dirac-picture. 

PACS numbers: 73.22.Pr, 71.70.Di, 81.05.uc, 71.70.-d 



I. INTRODUCTION 

Graphene nanostructuresi"— attract increasing at- 
tention mainly due to their potential applications in 
high mobility electronicsi^^ and sohd state quantum 
information processings^ In particular, low nuclear 
spin concentrations expected in graphene promise long 
spin lifetimesi^"— and make graphene quantum dots 
(QDs)i"— interesting for spin-qubit operations.— More- 
over, graphene nanostructures may allow to investigate 
phenomena related to massless Dirac Fermions in re- 
duced dimensions liii^r^ Intensive research has been trig- 
gered by the unique electronic properties of graphene^^ 
including the gapless linear dispersion, and the Lan- 
dau level (LL) spectrumi^"— Recent advances in fab- 
ricating width-modulated graphene nanoribbons have 
helped to overcome intrinsic difficulties in creating tun- 
neling barriers and confining electrons in graphene, 
where transport is dominated by Klein tunneling-related 
phenomenai^'^ Graphene QDs have been fabricated and 
Coulomb blockadejii^ quantum confinement^ and charge 
detection^ have been demonstrated. 




In this article, we focus on the eigenenergies of 
graphene quantum dots (see Fig. [1} as a function of a 
perpendicular magnetic field. In graphene, the linear 
band crossing at the so-called Dirac point suggests a 
close connection between the dynamics of electrons and 
free, ultrarelativistic Dirac particlesi^ One might there- 
fore expect a magnetic-field dependence of quantum dot 
eigenenergies that closely mirrors that of massless Dirac 
particles. Indeed, this connection has been used recently 
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FIG. 1: Shapes and sizes (50x50 nm) of graphene quantum 
dots confined by (a) a smooth valleyspin-conserving potential 
[Eq. (|lip . the length scale of the confinement is marked by 
Ad], (b) atomically sharp armchair and zigzag boundaries. 
Dots with disorder due to (c) bulk defects or (d) edge rough- 
ness. 



to discuss the spectrum of ideal, circular graphene dots 
with smooth confinement , ^^i^^ However, in more realistic 
models of finite graphene nanostructures, quantum con- 
finement, edge effects, and lattice defects introduce a host 
of competing length scales absent from the simple Dirac 
picture. Much progress has been made in understand- 
ing the unique LL spectrum, and the resulting Hall ef- 
fect, in graphenci^"— The magnetic-field dependence of 
the addition spectrum has been exploited in recent work 
to (approximately) pin down the electron-hole crossover 
point <^ In the present paper we report on a systematic 
study of the B-field dependence of electronic eigenstates 
of graphene quantum dots of experimentally realizable 
size (diameter d < 80 nm). We highlight the interplay 
of different length scales controlling the break-down of 
the valley symmetry by K — K' scattering. The latter 
is found to be key to the understanding of the diamag- 
netic spectrum. We find the i3-field dependence of the 
level variance to be a sensitive measure for the strength 
of K — K' scattering and obtain good agreement with 
experimental Coulomb blockade data. 

The paper is organized as follows: we first briefiy sum- 
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marizc the Dirac picture of Landau level formation for 
massless charged Dirac particles and discuss the length 
scales relevant to its applicability to finite-size graphcne 
quantum dots (Sec. II). In Sec. Ill we present realistic 
simulations for graphene quantum dots with zigzag and 
armchair edges, with edge roughness as well as with bulk 
disorder. A comparison between the calculated B-ficld 
dependence of the level variance and experimental data 
is given in Sec. IV, followed by a short summary (Sec. V). 



II. THE DIRAC PICTURE AND ITS 
LIMITATIONS 

The remarkable similarity of the low-energy band 
structure of graphene with the dispersion relation of 
a massless Dirac particle in two dimensions has been 
widely exploited in a variety of theoretical models for 
graphene.— However, the applicability of such models 
requires careful consideration of competing effects that 
go beyond the simple, yet intriguing Dirac picture42ii^ 
A case in point is the diamagnetism, i.e., the magnetic 
response of a finite-size graphene quantum dot. It is of 
considerable interest to inquire into the applicability as 
well as the limitations of the well-known diamagnetic the- 
ory of charged massless Dirac fermions. 

The magnetic-field {B) dependence of the spectrum of 
free Dirac particles was first solved in an early paper by 
Rabi^ shortly after the Dirac equation was proposed. 
The Dirac equation for a massless particle with charge 
q{= — |e|) in the presence of a potential V{x) with time- 
like coupling and a perpendicular, homogeneous mag- 
netic field B = V X A = —ByV x ex reads 

HD = Ho + HB=VF<j-[p-^A^+aoVix), (1) 

with a ~ (ax, ctj,) the Pauli matrices and do = 1- In the 
limiting case of strong magnetic field where IqA/cj 3> 
|y(x)| the solution of Eq. ([T]) predicts the formation of 
Landau levels A— 

E^{B) = sgn{n)^2\e\hvl\n\B, n G Zq. (2) 

We explicitly label this reference spectrum with the 
superscript "D" (for Dirac equation). Equation ([2]) 
contains several remarkable features absent from non- 
relativistic diamagnetism: a ground state Landau level 
n = the energy of which does not depend on B at 
all. Higher Landau levels n = ±1,±2,..., arc dis- 
tributed symmetrically around n = 0, and feature a \/B 
rather than a linear dependence on B known from non- 
relativistic diamagnetism. The high-field regime [Ec^. ([2])] 
is controlled by just two length scales, the (energy depen- 
dent) de Broglie wavelength Ap and the magnetic length 
Ib = y/hc/ (eB). The strong (weak) field regime is char- 
acterized by ^ Ap (^B ^ Ap). In the limit of weak 
magnetic fields, Eq. ([T]) predicts the lowest-order energy 



corrections to scale linearly with B, unlike the conven- 
tional non-relativistic behavior (cx B^). Eigenstates of 
Eq. ([1]) form two-spinors with definite helicity (or "chi- 
rality") 

M^> = ^'?-|jlV'>=±^l^) (3) 

in the absence of external fields. 

The ideal, infinitely extended graphene sheet featur- 
ing a honeycomb lattice made up by two interleafed 
triangular sublattices (A and B), can be described in 
nearest neighbor tight-binding approximation by the 
Hamiltonian^ 

where the sum extends over pairs of adjacent lattice 
sites, |(/)j,s) is the tight-binding orbital at lattice site j, 
Vi is a locally varying potential, and t (of the order of 
2.8 eV) is the nearest neighbor hopping matrix element. 
[In the numerical calculations we take into account sec- 
ond and third nearest neighbor coupling^M^iiS in addi- 
tion to Eq. Q in order to quantitatively account for the 
realistic band structure.] Close to the Fermi energy, the 
band structure of Eq. (|4|) can be approximated (assum- 
ing that Vi t) hy Si conical dispersion relation around 
the K point 

Eik + kK) = E{kK) + kdkE{kK) + O(fc^) « VF\k\, (5) 

where we have set E{kK) = 0. Note that the above ex- 
pansion ignores both the length scale of the graphene 
lattice constant a = 1.4 A and preferred directions of the 
lattice: due to the discrete lattice symmetry, the cone 
structure becomes squeezed along the K—K' directions, 
an effect known as triangular warpingi^i^ More impor- 
tantly, the band structure features two non-cquivalcnt 
cones ( "valleys" ) at the K and K' points in the reciprocal 
lattice. This additional degeneracy allows to formaly rep- 
resent the low-energy band-structure near iJ = in terms 

of Dirac-like four-spinors = (i/'a^^ ; V'b^^ > V'a^ ^ > ^) 
with amplitudes for the A-B sublattice in real space and 
for the K-K' valleys in reciprocal space. Operators in the 
four-spinor space can be represented by tensor products 
of ((To , (?) matrices acting on A-B sublattice amplitudes 
and analogous (tq , r ) Pauli-matrices acting on K ~ K' 
amplitudes. Choosing the origin in fc-space such that the 
connecting line between K and K' is along y, the effec- 
tive Dirac Hamiltonian in the absence of external scalar 
potentials becomes^ 

Hq = a ■ {p--A)®Ti+ a* ■ {p--A)®T2, (6) 
c c 

where ri^2 = {tq ± '^z)l'^- In addition to chirality, the 
valley-pseudospin projection 

r.|j)=j|j), J=±i (7) 
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associated with the vahey degree of freedom is conserved. 
The upper ( "particle- hke" , E > 0) and lower ("hole- 
like", E < 0) cones touching each other at K and K' 
with E = are related to each other by a particle-hole 
transformation 

C = cr^ ® To, CHC-^ = -H. (8) 

In the presence of a time-like scalar potential V{x)aQ (E) 
To, the Hamiltonian is invariant under an anti- unitary 
transformation ("time reversal"), T = iay (8) to C, where 
C denotes complex conjugation.— The wavefunctions at 
K and K' are related by time reversal symmetry. This 
symplcctic symmetry (T^ = — 1) is broken in the presence 
of a magnetic field 

fH{A)f-^ = H{~A) (9) 

lifting the two-fold Kramers-like degeneracy (Note that 
physical spin is not included in the present analysis.). 

We now consider a finite-size system of linear dimen- 
sion d, where l^(x) takes on the role of a confinement 
potential. With this additional length scale present, the 
Landau-level solution [Eq. is only valid in the strong 
magnetic field regime with Is d, while in the weak 
field regime. Is 3> d, the spectrum will be determined by 
V{x.). For zero magnetic field, eigenstates IV'if) and \ipK'} 
localized at the K and K' points are degenerate. Turning 
on a magnetic field lifts this degeneracy without (to low- 
est order) introducing couplings between K and K' . Fol- 
lowing first-order degenerate Rayleigh-Schrodinger per- 
turbation theory, the perturbation matrix W describing 
the lowest order correction to the field-free spectrum of 
Eq. ([6]) takes the form [for the gauge A = {—By, 0, 0)] 

W = {Hb)(Jcc®t,, (10) 

where {E[b)kk ^9 ("^sl 2/ IV'a) /c is real, and linear 
in B. The perturbation preserves the valley symmetry. 
The absence oi K — K' coupling and the linear magnetic 
field dependence of {Hb)xx ^'^'^ ca,c\\ decoupled Dirac 
cone implies that energy eigenvalues linearly cross the 
line i? = in pairs of two, forming an x-shaped intersec- 
tion (see Fig. [5]) . This is in contrast to the nonrelativistic 
diamagnetic response oc B^ in the perturbative limit. It 
rather resembles the paramagnetic level splitting in con- 
ventional quantum dots when the magnetic field lifts a 
degeneracy. Examples of the latter are lifting of Kramers' 
degeneracy^ or a symmetry-induced degeneracy as in 
circular quantum dotsj^ 

In order to quantitatively simulate the diamagnetic re- 
sponse of a finite-size graphene quantum dot we first con- 
sider a smooth confinement potential that is slowly vary- 
ing on a length scale of the lattice constant such as to 
approximately conserve the valley-pseudospin projection 
Tz (or valley symmetry). Moreover, it conserves particle- 
hole (anti-) symmetry [Eq. ([5])] in order to prevent Klein 
tunneling. Such a potential was first proposed by Berry 
and Mondragoni^, in the context of neutrino billiards, 

y(x) = K,(e'^'-(^)/^'^-l)a,, (11) 
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FIG. 2: (color online) (a) Magnetic-field dependence of 
the eigenenergies of a graphene quantum dot with smooth 
confinement which approximately preserves valley symmetry 
[Eq. (fTTj) . see Fig.flja)]. Landau levels for n = ±1, ±2 [dashed 
lines, see Eq. ((2)1] are inserted as guide to the eye. The four 
symbols (A, T, ■, ♦) mark parameter values for which eigen- 
states are shown in Fig.|3l (b) Close-up of the avoided crossing 
of two eigenstates in (a). Dots represent numerical data, the 
continuous line is a fit to Eq. (|10p . (c) Close-up of avoided 
crossings around the diabatic ridge formed by the first Landau 
level (see text). 



where Ar(x) is the outward distance from the quantum 
dot boundary and Ad introduces an additional length 
scale controlling the preservation of valley symmetry. 
We choose Ad ~ 24A [sec Fig. [IJa)] much larger than 
the lattice spacing (Ad ^ a sa 1.42A). Consequently, 
Eq. ([TT]) varies slowly on the scale of the lattice con- 
stant, conserves valley symmetry [to order (a/Ad)^] and 
provides a realization of the (approximately) K — K' de- 
coupled diamagnetic perturbation [Eq. (fTOj) ]. We note 
that realizations of potentials of the form of Eq. ([TT]) 
are, to our knowledge, currently experimentally not avail- 
able. We employ a third-nearest neighbori^SiiSid^ tight- 
binding approximation (to correctly describe triangular 
warping) and simulate a 50x50 nm graphene QD con- 
taining « 100.000 carbon atoms. The magnetic field is 
included by a Peierls phase factor. We use a Lanczos 
diagonalization in conjunction with an LU factorization 
to efficiently calculate the 500 eigenvalues closest to the 
Fermi edgo^— [see Fig. [5]. 

In the limit of weak magnetic fields, we find that our 
numerical results, indeed, follow the linear _B-field depen- 
dence of the energy eigenvalues as predicted by pcrturba- 
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FIG. 3: (color online) Typical eigenstates (plotted is the ab- 
solute square of the wavefunction) of a graphene quantum dot 
with smooth confinement [see Fig. [1] (a) and Eq. at high 
magnetic field (B=25 T), corresponding to the zeroth [(a) and 
(b)] and first [(c) and (d)] Landau level. Symbols (A, ▼, ■, ♦) 
correspond to those marking the position of these states in 
the energy level diagram [Fig. [2ja)]. 



tion theory [Eq. ([T0| ] [see Fig.[2l^b)]. Residual deviations 
from the perfect lattiee symmetry (due to the finite width 
Ad of the confinement) and, thus, weak non-conservation 
of Tz appear as minute energy splittings between near- 
degenerate levels in the B ^ limit. The resulting level 
splitting at B = is, however, very small (120/i eV), i.e., 
two orders of magnitude below the mean level spacing 
(« 10 meV). 

Turning now to the high-field limit, Ib ^ d, the influ- 
ence of conflnement effects should be diminished and the 
formation of Landau levels following the Dirac picture 
[Eq. ([!])] is expected. The transition from low to high 
magnetic flelds drastically changes the density of states 
(DOS). The depletion of the DOS near £; = at low 
flelds, 



p{E) 



2{hvFy 



\E\ 



(12) 



is replaced, for increasing B, by an increasing number 
of eigenstates moving towards the Landau level at Eq , 
which is located at the Dirac point [see Fig. [JJa)]. More 
speciflcally, all graphene levels with energies in between 
the two first Landau levels, E^i < E < ijf , adiabati- 
cally converge to the level at Eq =0. As we have shown 
recently,— this unique feature can be used to pin down 
the energetic position of the Dirac point in the exper- 
iment and thus of the electron-hole crossover region in 
real graphene quantum dot devices. 

While, at low fields, the valley symmetry is approx- 
imately preserved by the potential in Eq. (jll[) . a large 



number of sizeable avoided crossings appear at higher 
magnetic fields as the magnetic length is reduced to 
Ib ^ d. Edge states that couple to bulk states or to other 
edge states become prevalent. The complicated pattern 
in Fig. [2ja),(c) of many avoided crossings near the first 
Landau energy E = E-^i reflects this interplay between 
magnetic bulk and edge states. Levels with eigenen- 
ergies that follow the predicted values for the Landau 
levels, E^, are localized in the interior of the quantum 
dot and well separated from the edges. Conversely, the 
states with energies in between the values E^ should 
be strongly influenced by the spatial conflnement in the 
quantum dot. Wave functions corresponding to energy 
levels close to n = and n = l Landau levels [Fig.[3l[a,c)] 
as well as those in between n = and n = ±1 [Fig. [3Jb)] 
and between n = 1 and n = 2 [Fig. EJd)] confirm these 
expectations. 

Apparently, the typical level splittings at the avoided 
crossings [Fig. [Hb,c)] are dramatically enhanced in the 
high-field regime. This is due to the fact that, as com- 
pared to the low-field case, the amplitudes of wave func- 
tions are enhanced at the dot boundary [see Fig. [31[b)]. 
Since, in addition, the lattice symmetry is broken at the 
boundary, these edge states do have an increased cou- 
pling strength to all other states in the sample. Fol- 
lowing the Wigner-von Neumann non-crossing rule^ this 
increased coupling strength leads to increased level split- 
tings at the crossing point. 



III. REALISTIC GRAPHENE QUANTUM DOTS 

A. Clean dots with sharp edges 

We turn now to realistic graphene quantum dots where 
the nanostructures are terminated by atomically sharp 
edges, either of armchair or zigzag shape [Fig.[TJb)]. Fol- 
lowing recent estimates for passivation of the dangling 
carbon bonds at the edges of graphene samples (e.g., by 
attached hydrogen)^ we set the potential of the out- 
most carbon atoms to 1.7 eV. Although in this case many 
of the surface states present in a perfect zizag bound- 
ary remain suppressed, such a confining potential leads 
to substantial changes in the energy level spectrum [see 
Fig.m^a)] as compared to the model potential [Eq. ([TT|) ]. 
Most importantly, in the low-field regime [see Fig. HJb)] 
the linear _B-field dependence is replaced by a quadratic 
dependence (cx _B^) of the level splitting resembling the 
nonrelativistic diamagnetic response. This is due to the 
presence of sizable avoided crossings near B = as a 
result of broken valley symmetry caused by the edges. 
The valley symmetry is, thus, broken upon refiection at 
atomically sharp "clean" zigzag edges and Tz is no longer 
conserved. In terms of perturbation theory, the confining 
potential V now includes off-diagonal components in the 
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FIG. 4: (color online) Same as Fig. [J] but for a quantum dot 
with atomically sharp zigzag and armchair boundaries. The 
solid line in (b) is a fit to Eq. (|14|l . the dashed line (corre- 
sponding to Vkk' ~ 0) is inserted as guide to the eye. The 
evolution of one eigenenergy with magnetic field is drawn with 
a thick green line as guide to the eye in (a), (c). The arrow 
in (c) marks a kink in the magnetic field dependence of this 
state (see text). 



following the first Landau level. We rather observe a bun- 
dle of states sequentially moving along the characteristic 
energy of the first Landau level, much like in a relay 
[see Fig. IH^c)] . Such an interrelated sequence of avoided 
crossings is well-known from atomic physics as "diabatic 
ridge" - riding states localized on potential barriers.— A 
direct consequence is that the evolution of eigcnstates for 
an increasing magnetic field [see green highlighted line in 
Fig. m^a)] features sharp "kinks" when crossing the ridge 
following the first Landau level [see arrow in Fig. |H[c)] . 
As the state is transiently trapped by the ridge, it moves 
away from the Dirac point, and continues again mono- 
tonically towards the Dirac point once clear of the ridge. 
These kinks due to the ridge riding mechanism have been 
observed in the experiment^ serving as an additional in- 
dicator for the position of the lowest Landau level and 
the clcctron-holc crossover. 

The present results show that the atomically sharp 
edges destroy the linear i3-dependence of the energy lev- 
els at weak fields, but they do preserve the square-root in- 
dependence at very high fields. The linear B dependence 
results from the fragile suppression of K — K' scattering 
(cx {V)j^j^,) between the Dirac cones while the square- 
root dependence results from the much more robust dis- 
persion relation of the individual cone. Therefore, Lan- 
dau levels survive the introduction of sharp boundaries 
much better than the energy levels at weak fields. In 
turn, even when in the experiment a Dirac-likc Landau 
level spectrum is recorded, many other features of the 
same graphene sample may very well show large devia- 
tions from predictions based on Dirac theory. 



pseudospin degree of freedom 

W ^C7^(E) {Hb)kk '^z + Re {V)j^j^, + Im (V")^^, r^, 

(13) 

with coupling matrix elements between the valleys 

{V)kk' = <^o® (V'U^IV')^' + <^oc ® (V'U^I'/'>i3' and 
eigenvalues 

e^±^{HB)lK + {^e/2)\ l:,e ^ 2\{Vkk')\ ■ (M) 

The coupling between K and K' cones thus leads, accord- 
ing to the Wigner-von Neumann non-crossing rulc)^ to 
avoided crossing with level splittings Ae [see Fig. HJ^b)] 
proportional to the coupling strength Vkk' between the 
two Dirac cones. Conversely, a fit to Eq. (|14p yields a 
sensitive indicator for the amount of K—K' scattering in 
the quantum dot. 

For high magnetic fields [sec Fig. |4l[c)], the presence 
of Vkk' coupling lead to a large number of correlated 
avoided crossings when the edge states move towards the 
zeroth bulk Landau level. In other words avoided cross- 
ings appear when the energy of eigenstates evolving to- 
wards the level "pass" through the energy E^{B) of 
the first Landau level [Fig. IH^c)]. Due to a large num- 
ber of avoided crossings, there are no states continuously 



B. Dots with bulk disorder 

To further elucidate the role of lattice symmetry break- 
ing in graphene quantum dots, we now consider isolated 
lattice defects in the bulk. The preceding results sug- 
gest that disorder realizations that break the K — K' 
symmetry arc of crucial importance. We therefore con- 
sider first single lattice vacancies [see inset in Fig. [SJa)] 
with defect densities of rii = 10~^ to 10~^ impurities 
per carbon atom. To isolate the effect of such bulk de- 
fects from the K — K' scattering of the edges, we first 
use the smooth symmetry preserving boundary poten- 
tial [see Eq. (|lip ]. Overall, the diamagnetic spectrum 
closely resembles that of clean samples with atomically 
sharp edges (see Fig. [5]). In particular, avoided crossings 
near B = with a quadratic field dependence as well as 
the formation of sequences of avoided crossings along the 
ridges of bulk Landau levels are found. Wave functions 
of eigenstates at these energies display patterns which 
are very similar to those for the clean system [compare 
Fig. ^h) and Fig. IH^a)]. Likewise, the kink pattern at 
the crossing of the ridges appears robust against disorder 
[sec, e.g.. Fig. [Sljc)]. We do, however, observe new ridges 
between the Landau level energies which were absent for 
edge scattering [see dotted lines in Fig.[51Ja,c)] . The corre- 
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FIG. 5: (color online) Same as Fig. [2] for a quantum dot fea- 
turing 30 single lattice vacancies (see inset) out of a total 
number of 100.000 carbon atoms. Dotted lines mark the evo- 
lution of two states localized at one defect. Symbols mark the 
states for which the corresponding wavefunction is shown in 
Fig.H 

spending eigenstates near these new ridges [see Fig. (SJ^b) , 
marked by ▼ in Fig. [S{a)] are pinned to a single defect, 
where the lattice periodicity and the A-B sublattice sym- 
metry are broken, and the partitioning of the wavefunc- 
tion in four components according to the Hamiltonian 
in Eq. ^ fails. Such localized states can therefore be 
expected to behave differently from the bulk Landau lev- 
els. The resulting ridges feature a very weak quasi-linear 
magnetic field dependence. We therefore conjecture that 
these structures are due to avoided crossings with such 
localized defect states. Recent analysis^^*^ has shown 
that K~K' scattering at impurities in graphene, i.e. the 
strength of Vkk' ^ tends to be strongly energy dependent. 
Bound states due to adsorbates (leading to an enhanced 
local density of states at the adsorption site) frequently 
have energies close to the Dirac point i^iS Since the lo- 
calized defect states with energy En feature a very weak 
explicit linear magnetic field dependence, 

mcV 

En[VKK',B]^En[VKK']+a\B\, |a| < 1^, (15) 

the implicit magnetic field dependence of VKK'{En[B]) 
may be neglected. While the detailed energy depen- 
dence of Vkk' may be connected to the specific defect 
present (e.g.. Stone- Wales defects, or attached nitrogen 
molecules) , we still qualitatively expect an analogous lin- 




FIG. 6: (color online) Eigenstates of a graphene quantum dot 
with disorder: (A, T) 30 single vacancies [one single vacancy 
shown as inset in Fig.[5{a)] or (■, ♦) edge roughness of ±2 nm. 
Positions are marked by corresponding symbols in Fig. [Sfa) 
and Fig. [Ha). 



ear i?-field dependence. Such localized states with weak 
magnetic field dependence were also observed experimen- 
tally in Coulomb blockade measurements 

In the low field regime, we observe a significant change 
in the B-evolution of eigenenergies: avoided crossings 
become asymmetric [see Fig. El^b)]. The reason is that 
the K — K' splitting introduced by the lattice vacancies 
is strong enough to yield different matrix elements for 
(-^)kk ^^^'^ (^)k'K'- consequently, the slope of both 
eigenvalues of Eq. (fT5|) is different [see Fig. [SJb)]. To il- 
lustrate that the diamagnetic spectrum, in particular the 
avoided crossing distribution, is due to the breaking of 
the A-B sublattice and, in turn, to the breaking of valley 
pseudospin symmetry induced by the defects, we present 
as counter example the spectrum for double vacancies. 
We introduce such a vacancy in accordance to our soft- 
wall potential, Eq. (|lip . As such double vacancies act on 
an entire unit cell in the hexagonal lattice, they approxi- 
mately conserve the electron-hole and K—K' symmetry. 
We find, indeed, that for the same number of defects as 
of single vacancies (Fig. [5]), only avoided crossings with 
comparatively small energy splittings appear in the en- 
ergy level diagram (Fig.[7|), resembling much more closely 
Fig. [3] This clearly indicates that it is the breaking of 
the A-B symmetry and not the presence of defects per 
se which is responsible for the break-down of the Dirac 
picture for graphene quantum dots. 
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FIG. 7: (color online) Same as Fig. [2] for a quantum dot 
featuring 30 double lattice vacancies (see inset) in 100.000 
carbon atoms. 



C. Dots with edge disorder 

We consider now a clean graphene quantum dot with 
atomically sharp but disordered edges. We connect short 
straight edge segments (with a random length between 
0.5-3 nm) to obtain a polygon-shaped boundary [see 
Fig. [TJc)] with a disorder amplitude of \Sd\ < 2 nm. 
6d is thus (at energies close to the Dirac point) smaller 
than the wavelength Ap of the confined particles as well 
as the magnetic length, but larger than the lattice con- 
stant. Since rough boundaries, just like bulk defects, lo- 
calize stateS)^i2i we expect similar signatures of these 
two types of disorder. Indeed, in the low field regime, 
dots with rough edges feature a similar pattern of fluc- 
tuating energy levels as dots with single-lattice defects 
(compare Fig. |S] and Fig. [5]). The spectra are so sim- 
ilar that it is difficult to distinguish between bulk and 
edge disorder breaking the K — K' symmetry. Also in 
the high-field regime, the evolution towards the Landau 
levels features the correlated sequence of avoided cross- 
ings reflecting the diabatic ridges [see Fig. [5Jc)]. Wave 
functions of cigenstatcs at these energies display patterns 
which are very similar to those for the clean system [com- 
pare Fig. [31[a) and Fig. IHl^d)]. In particular, states near 
Landau levels localize in the interior of the dot, and thus 
are not influenced by edge disorder. Likewise, the kink 
pattern at the crossing of the ridges is robust against edge 
disorder. While the inclusion of edge disorder does not 
give rise to qualitatively new effects on the eigencnergy 
spectrum in the high-field regime, differentiating between 
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FIG. 8: (color online) Same as Fig. |4] for a quantum dot fea- 
turing rough edges (with a roughness amplitude Sd = ±2 nm). 
Symbols mark the states for which the corresponding wave- 
function is shown in Fig. [S] 
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FIG. 9: (color online) (Avoided) crossings for (a) soft edges, 
(b) hard edges, (c) rough edges, and (d) rough edges plus bulk 
disorder. The level splitting is (a) 0.1 meV, (b) 2 meV,(c) 
3 meV, and (d) 2.5 meV. 



edge and bulk disorder might become possible by probing 
the different scaling behavior for bulk and edge disorder 
with the size of the graphene quantum dot. 

It is instructive to directly contrast the level splitting 
Ae of an avoided crossing [Eq. due to finite K~K' 
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coupling for different scenarios: (i) valley-symmetry pre- 
serving confinement [Fig.[3l^a)], (ii) clean graphene quan- 
tum dots with atomically sharp edges [Fig. EJb)] and 
(in) disordered graphene quantum dots [Figs. ElJCjd)]. 
We observe an increase in the size of the average level 
splitting (Ae) due to K — K' coupling at the edges and 
impurities. In scenario (i) (Ae) is at least one order 
of magnitude smaller than the mean level spacing ((5e), 
(Ae) ^ {5e). In (ii) where (Ae) < (fe) a first-order per- 
turbative treatment of K—K' coupling correctly describes 
the level repulsion. For localization at defects (iii) where 
(Ae) < {5e) level pairs at K and K' feature different 
quadratic dependences on B, i.e., they can no longer be 
parametrized by Eqs. (fT3|) and (fT4|). 



IV. 



(a) 2. 



COMPARISON WITH EXPERIMENT: 
LEVEL SPACING FLUCTUATIONS 



For a comparison with the experimental data for the 
magnetic response of graphene quantum dots we pur- 
sue two strategies: in a direct approach we compare our 
models with the observed parametric B-field evolution 
of individual Coulomb blockade peaks. Alternatively, we 
identify the B-field dependence of the level spacing fluc- 
tuations (variance) as a robust measure for the degree 
of disorder in a graphene quantum dot. Specifically, we 
determine the rescaled (or unfolded) variance (T^ of the 
distribution of neighboring energy level spacings 5e 



(fc) 



(16) 



Since switching on a magnetic field B leaves the number 
of states unchanged, {5e) is (approximately) independent 
of _B, while higher moments of the level distribution are 
drastically affected. At i? = 0, pairs of energy levels 
are split by the characteristic energy Ae of the avoided 
level crossings. However, as long as the mean width of 
the avoided crossing (Ae) [Eq. p^ ] is small compared to 
the mean level spacing {5e) , the level sequence fluctuates 
between small and large spacings, while spacings of the 
order of Se are unlikely. We thus expect for magnetic field 
B — Q a. comparatively large variance a^- For increasing 
\B\ the levels become more equally spaced, leading to 
a decrease in ae- Correspondingly, the variance a. of 
the level spacings should feature a peak at i? = 0. The 
numerical results for the dependence of on B for dif- 
ferent disorder strength (i.e., different [Vftrifl) are shown 
in Fig. 1101 Our data display a pronounced peak of 
for the clean flake slowly decreasing for increasing num- 
ber of single- vacancy defects (i.e., for increasing K — K' 
scattering) [see Fig. [TUT a)]. Note that both the peak 
height bX B = Q and the overall value of decrease 
with increasing disorder. The latter can be explained 
by the emerging localized states that feature a regular 
spacing (and hence a suppressed variance CTe). Conse- 
quently, if a given spectral region is more prone to fea- 
ture localized states due to adsorbates than others i^^i^^ 






B[T] 



FIG. 10: (color online) Variance of the mean level spacing 
[see Eq. {TS])] as a function of magnetic field for (a) single- 
vacancy defects and (b) double-vacancy defects (note the ex- 
panded y-scale) for three values of disorder concentrations: 
from top to bottom = 3 ■ 10"^ 6 • 10"* and 2 ■ 10"*. 



we expect an energy-dependence on a depending on the 
specific type of adsorbate. Such an energy dependence 
could be exploited to directly measure the energy depen- 
dence of X — /^'scattering by determining the statistics 
of Coulomb blockade resonances at different back-gate 
voltages. For comparison we also show for double va- 
cancies preserving A-B symmetry, [see Fig. Iior b*). note 
the different scales]. Accordingly, cr^ is, indeed, strongly 
dependent on the amount of K — K' scatterers, not on 
the overall number of defects. 

The decrease in peak height with increasing K — K' 
scattering should thus provide a robust and sensitive 
measure for K—K' scattering present in the experiment. 
To test this conjecture, we have measured the evolution of 
42 Coulomb blockade peaks for varying magnetic fieldi^ 
We follow the parametric motion of the peak positions, 
and hence the eigenstate energies, as a function of mag- 
netic field B € [—2, 2]T. To compare with our numerical 
results, we take into account a charging energy of 13 meV 
(determined independently^) as well as spin (by Zeeman 
splitting). We observe, indeed, a quadratic i3-field de- 
pendence rather than the linear dependence predicted for 
conserved K — K' symmetry of the Dirac equation. Our 
experimental data can be well described by Eq. [see 
Fig. [TlTa')]. For pairs of consecutive Coulomb blockade 
resonances belonging to the same avoided crossing, we 
find a mismatch in slopes, in agreement with our numer- 
ical findings for the rough-edged quantum dot [compare 
Fig. fTlTa) and Fig. in^c,d)]. This has important conse- 
quences for the interpretation of experimental data: the 
roughness present in the experimental dot does not allow 
to disentangle K and K' states. To be more quantita- 
tive, we compare our simulations for the level variance 
with experimental data. We indeed find good agreement 
as confirmed by a noticeable peak also in the experimen- 
tal data for erg (B) [see Fig. lTlT b)]. The offset between the 
two data sets in Fig. [TUb) is attributed to a possible en- 
ergy dependence oiK—K' scattering as well as statistical 
fluctuations in charging energy and in the number of lo- 
calized states for different values of the back gate voltage. 
By using the edge roughness 5d as the only adjustable 
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FIG. 11: (color online) (a) Coulomb blockade peaks measured 
as a function of applied plunger gate electrode potential (Vg) 
and mag netic field applied perpendicular to the sample (see 
Ref. [50] for details about the measurement). Lines: fit to 
Eq. p4|) for two measured Coulomb blockade peak pairs with 
a level splitting at B = of ~ O.SV. (b) Normalized variance 
of the level spacing o-g [see Eq. (|16p ] as a function of magnetic 
field. Dots denote an average over 20 consecutive experimen- 
tal Coulomb blockade peaks for two values of back gate volt- 
age Vbg [corresponding to the upper {Vbg ~ 20 V) and lower 
{Vbg ~ 38 V) data points]. Solid lines represent simulations 
for different edge roughness amplitude 5d [Fig. [Ud)] as only 
fit parameter: 5d ~0.5 (0.6) nm, for the upper (lower) curve. 



parameter, we can match the measured (Je{B) very well 
with our numerical simulation [see lines in Fig. Illf b)] . 
Good agreement is found for an edge roughness of about 
6d « 0.5 ± 0.2 nm, or, equivalently, a K — K' scatterer 
concentration in the bulk rii « (3.5 ± 1) • 10"'', both of 
which are well within expectation. We emphasize that, 
although we can quantify the resulting overall strength of 
K—K' coupling in our experimental quantum dot, we can- 
not disentangle whether the observed K — K' scattering 
comes from edge roughness, lattice defects, or disorder 
through flake-substrate interactions with a length scale 
comparable to the lattice constant. 



sity of states to the emergence of Landau levels. At 
high field strength, we find that Landau levels follow 
the square-root dependence of the Dirac equation, mani- 
fested in the energy level diagram by sequences of corre- 
lated avoided crossings along "diabatic ridges". These 
ridges lead to characteristic kinks in the evolution of 
energy states that cross a Landau level. Appearing 
also in Coulomb blockade measurements, these kinks can 
be used to experimentally pin down the electron-hole 
crossover point In the perturbative regime of small 
magnetic fields, we find that the linear dependence on B 
predicted by the model of massless Dirac fermions disap- 
pears when the valley symmetry is broken. Even perfect 
armchair and zigzag edges are suflacient to break the sub- 
lattice symmetry giving rise to avoided crossings with a 
quadratic dependence on B instead. A similar effect is 
observed for lattice defects: single lattice vacancies break 
the K — K' symmetry and thus result in avoided cross- 
ings with substantial level splittings (even for defect con- 
centrations as low as 1 in 20.000). By comparison with 
double vacancies which conserve the sublattice symme- 
try we show that it is not the presence of disorder per se 
which leads to deviations from predictions by the Dirac 
equation, but the breaking of valley symmetry. 

We compare our theoretical predictions with exper- 
imental results on the parametric B-field evolution of 
Coulomb blockade peaks. As a quantitative indicator 
for the strength of K — K' scattering, we identify the 
variance of the level spacing distribution. We observe 
a peak in the variance at B = due to level correla- 
tions near avoided crossings. We find quantitative agree- 
ment between the measured and the calculated data for 
CTg which enables us to pinpoint the amount of K — K' 
scattering present in our experimental flake. The present 
results provide a sensitive indicator for the quality of the 
graphene dot and demonstrate the limits of the Dirac 
picture in describing the experiment. 
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